The Baikal Deep Underwater Neutrino Experiment: 
Results, Status, Future * 
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We review the present status of the Baikal Underwater Neutrino Experiment and present results obtained with the 
various stages of the stepwise increasing detector: NT-36 (1993-95), NT -12, (1995-96) and NT-96 (1996-97). Results 
cover atmospheric muons, first clear neutrino events, search for neutrinos from WIMP annihilation in the center of 
the Earth, search for magnetic monopoles, and - far from astroparticle physics - limnology. 
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The Baikal Neutrino Telescope is being deployed in Lake Baikal, Siberia, 3.6 km from shore at a 
depth of 1.1 km (see fig.l). At this depth, the light absorbtion length for wavelengths between 470 
and 500 nm is about 20 m. Scattering is strongly forward peaked ((cos 9)) 0.95). Typical values 
for the scattering length are about 15 m. Expressed in terms of the effective scattering length 
L e ff = L scat t/(l — (cos 6)), this corresponds to L e ff = 300 m. 



NT-200, the medium-term goal of the collaboration |y, |2j, will be finished in April 1998 and will 
consist of 192 optical modules (OMs). An umbrella-like frame carries 8 strings, each with 24 
pairwise arranged OMs. Three underwater electrical cables connect the detector with the shore 
station. Deployment of all detector components is carried out during 5-7 weeks in late winter when 
the lake is covered by thick ice. 

In April 1993, the first part of NT-200, the detector NT-36 with 36 OMs at 3 short strings, was put 
into operation and took data up to March 1995. A 72-OM array, NT- 72, run in 1995-96. In 1996 
it was replaced by the four-string array NT-96. Summed over 700 days effective life time, 3.2 • 10 8 
muon events have been collected with NT-36, -72, -96. Since April 6, 1997, NT-144, a six-string 
array with 144 OMs, is taking data (see fig. 2). 

The OMs are grouped in pairs along the strings. They contain 37-cm diameter QUASAR PMTs. 
The two PMTs of a pair are switched in coincidence in order to suppress background from biolu- 
minescence and PMT noise. A pair defines a channel. 



* Talk given at the Int. School of Nuclear Physics, Erice, Sept. 1997 
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Figure 1: The Baikal detector complex (status since 1994). 1,2 - wire cables to shore, 3 - opto-electrical 
cable to shore, 4,5,6 - string stations for shore cables 1,2,3, respectively, 7- string with the telescope, 8 - 
hydrometric string, 9-14 - ultrasonic emitters 




Figure 2: Schematic view of the Baikal Telescope NT-200. The modules of NT-144, operating since April 
1997, are in black. The expansion left-hand shows 2 pairs of optical modules ("svjaska") with the svjaska 
electronics module, which houses parts of the read-out and control electronics. Top right the array NT-96 
is shown, which took data between April 1996 and March 1997. 
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A muon-trigger is formed by the requirement of > N hits (with hit referring to a channel) within 
500 ns. N is typically set to the value 3 or 4. For such events, amplitude and time of all fired 
channels are digitized and sent to shore. The event record includes all hits within a time window 
of -1.0 usee to +0.8 //sec with respect to the muon trigger signal. A separate monopole trigger 
system searches for clusters of sequential hits in individual channels which are characteristic for the 
passage of slowly moving, bright objects like GUT monopoles. 

In the initial project of NT-200, the optical modules were directed alternately upward and down- 
ward, with a distance of 7.5 m between pairs looking face to face, and of 5 m between pairs arranged 
back to back. Due to sedimentation of biomatter deterioiating the sensitivity of upward looking 
OMs we were forced to direct the OMs of the present arrays essentially downward. 

2 Track Reconstruction 

The parameters of a muon track crossing the detector are determined by minimizing § 



Here, ti are the measured times and Tj the times expected for a given set of track parameters and 
simplifying the theoretical picture to that of a "naked" muon not accompanied by electromagnetic 
showers. N^u is the number of hit channels, an are the timing errors. A set of parameters defining 
a straight track is given by 9 and (ft - zenith and azimuth angles of the track, respectively, uq and 
vq - the two coordinates of the track point closest to the center of the detector, and to ~ the time 
the muon passes this point. 

Only events fulfilling the condition "> 6 hits at > 3 strings" (trigger 6/3) are selected for the 
standard track reconstruction procedure which consists of the following steps: 

1. A preliminary analysis including several causality criteria which reject events violating the 
model of a naked muon. After that, a 0-th approximation of 9 and <j> is performed. 

2. A minimum search (minimization of based on the model of a naked muon. 

3. Quality criteria to reject most badly reconstructed events. 

The reconstruction procedure is described in detail in jj], [|]. Fig. 4 shows a typical single muon 
event firing 7 of the 18 channels of NT-36 and reconstructed with a x 2 /NDF = 0.57 (left), as well 
as the muon intensity 1^(9) as a function of measured zenith angle at the detector depth (right). 
The angular distribution is well described by MC expectations. 

3 Separation of Neutrino Events 

The canonical signature of neutrino induced events is a muon crossing the detector from below. 
With the flux of downward muons exceeding that of upward muons from atmospheric neutrino 
interactions by about 6 orders of magnitude, a careful reconstruction is of prime importance. Clear 
neutrino signals in the rather small NT-36 could be separated only over a very limited cone around 
the opposite zenith. Contrary to that, NT-96 can be considered as a real neutrino telescope for a 
wide region in zenith angle 9. 
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Figure 3: Left: Single muon event recorded with NT-36. Hit channels are in black. The thick line gives 
the reconstructed path, thin lines pointing to the channels mark the path of Cherenkov photons as given by 
the fit to the measured times. The sizes of the ellipses are proportional to the recorded amplitudes. Right: 
Zenith angle distribution of the muon intensity 1^(9) at a depth of 1170 m. Full circles: experimental data, 
solid line: 1^(9) calculated for stochastic energy loss. 

3.1 Search for Nearly Upward Moving Neutrinos with NT-36 

For separation of nearly vertical upward muons in the NT-36 prototype array we used special 
separation criteria instead of full reconstruction (see also ||). These criteria make use of two facts: 
firstly, that the muons searched for have the same vertical direction like the string; secondly, that 
low-energy muons generate mainly direct Cherenkov light and, consequently, are not visible over 
large distances and should produce a clear time and amplitude pattern in the detector - mostly 
only along one string. We have chosen the following criteria: 

1. \(U — tj) — (Tj — Tj) | < dt: the time pattern must be close to that of a straight 
upward moving muon. U(tj) are the measured times in any hit channels Ti(Tj) are 
the "theoretical" times expected for naked, up-going vertical muons and dt is a time cut. 

2. dAij(down — up) > 0.3: All down-looking channels must see clearly more light 
than any upward looking channel. dAij (down — up) = (Ai(down) — Aj(up)) / (Ai(down) + 
Aj(up)) and Ai(down)(Aj(up)) are the amplitudes of channel i(j) facing downward (upward). 

3. Ai(down)>Ape: Amplitudes of downward channels must exceed 4 photoelectrons 

4. dA(down — down) < 0. The light must not decrease from bottom to top, as one 
expects for showers generated by downward muons below the array. dA(down — 
down) is defined as that of the 3 possible combinations dA+j (down — down) = (Aj — Aj)/ (Ai + 
Aj) |j>j of downward channels which has the largest absolute value. For events due to showers 
below the array it peaks at values close to 1, for vertical neutrino candidates it should be 
close to zero. The criterium rejects nearly all shower events but only half of the neutrino 
sample. 

We analyzed the data taken with NT-36 between April 8, 1994 and March 5, 1995 (212 days 
lifetime). There were 6 PMT pairs along each of the 3 strings. The orientation of the channels 
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from top to bottom at each string was up- down-up- down-up- down. Upward-going muon candidates 
were selected from a total of 8.33 • 10 7 events recorded by the muon-trigger "> 3 hit channels". 
The samples fulfilling trigger conditions 1, 1-2, 1-3 and 1-4 with time cut dt = 20ns contain 131, 
17 and 2 events, respectively. Only two events fulfill trigger conditions 1—3 and 1~4- The first is 
consistent with a nearly vertical upward going muon and the second with an upward going muon 
with zenith angle 9^ = 15° (fig. 4). A detailed analysis || yields a fake probability of <3% for both 
events. 




Figure 4: The two neutrino candidates in NT-36. The hit PMT pairs (channels) are marked in black. 
Numbers give the measured amplitudes (in photoelectrons) and times with respect to the first hit channel. 
Times in brackets are those expected for a vertical going upward muon (left) and an upward muon passing 
the string under 15° (right). 

Considering the two neutrino candidates as atmospheric neutrino events, a 90 % CL upper limit of 
1.3 • 10 -13 (muons cm -2 sec" 1 ) in a cone with 15 degree half-aperture around the opposite zenith 
is obtained for upward going muons generated by neutrinos due to neutralino annihilation in the 
center of the Earth. The limit corresponds to muons with energies greater than the threshold 
energy E t ^ ~ 6 GeV, defined by 30m string length. This is still an order of magnitude higher than 
the limits obtained by Kamiokande Q, Baksan and MACRO Q. The effective area of NT-36 for 
nearly vertical upward going muons fulfilling our separation criteria 1-3 with cft=20ns is S e ff = 50 
m 2 /string. A rough estimate of the effective area of the full-scale NT-200 (with eight strings twice 
as long as those of NT-36) with respect to nearly vertical upward muons gives S e ff ~ 400 — 800 
m 2 . 

3.2 Fully reconstructed neutrinos separated with NT-96 

The NT-96 data are analyzed using the standard reconstruction procedure described in sect. 2 as 
well as a method similar to that described in the previous subsection. In this subsection we present 
results obtained with the standard procedure (see also Q). For NT-96, the most effective quality 
cuts are the traditional x 2 cut, cuts on the probability of non-fired channels not to be hit, and 
fired channels to be hit (P no hit and Phit> respectively), cuts on the correlation function of measured 
amplitudes to the amplitudes expected for the reconstructed tracks, and a cut on the amplitude 
X 2 defined similar to the time \ 2 defined above. To guarantee a minimum lever arm for track 
fitting, we reject events with a projection of the most distant channels on the track (Zdist) below 35 
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meters. Due to the small transversal dimensions of NT- 96, this cut excludes zenith angles close to 
the horizon, i.e., the effective area of the detector with respect to atmospheric neutrinos is decreased 
considerably (fig. 5, left). 

The efficiency of all criteria has been tested using MC generated atmospheric muons and upward 
muons due to atmospheric neutrinos. 1.8 • 10 6 events from atmospheric muon events (trigger 6/3) 
have been simulated, with only 2 of them passing all cuts and being reconstructed as upward going 
muons. This corresponds to S/N ~ 1. Rejecting all events with less than 9 hits, no MC fake event 
is left, with only a small decrease in neutrino sensitivity. This corresponds to S/N > 1 and the 
lowest curve in fig. 5. 

With this procedure, we have reconstructed 5.3 • 10 6 events taken with NT-96 in April/May 1996. 
The resulting angular distribution is presented in fig. 5 (right). 




Figure 5: Left: Effective area for upward muons satisfying trigger 9/3; solid line - no quality cuts; dashed line 
- final quality cuts; dotted line - final quality cuts and restriction on Z^ist (see text). Right: Experimental 
angular distribution of events satisfying trigger 9/3, all final quality cuts and the limit on Z^ist (see text). 



From the time period between April 16 and May 17, 1996 (18 days lifetime), three neutrino can- 
didates have been separated (see fig. 5, right), in good agreement with the expected number of 
approximately 2.3. Fig. 6 displays one of the neutrino candidates. Top right the times of the hit 
channels are shown as a function of the vertical position of the channel. At each string we observe 
the time dependence characteristically for upward moving particles. The angle regions ip min — ip max 
consistent with the observed time differences Aijj between two channels i, j are given by 

r ■ At - ■ 

cos(V> mm + v)< cos ^- < zo^ max - tj) (2) 

Tj -n 

with fi,fj being the coordinates of the channels, ip the muon angle with respect to fj — fi and r\ 
the Cherenkov angle. The bottom right picture of fig. 6 shows that the overlap region of all channel 
combinations of this event clearly lay below horizon. 

The same holds for the other two events, one of which is shown in fig. 7a. Fig. 7b, in contrast, shows 
an ambiguous event giving, apart from the upward solution, also a downward solution. This event 
is assigned to the downward sample. 
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Figure 6: A "gold plated" 19- hit neutrino event. Left: Event display. Hit channels are in black. The thick 
line gives the reconstructed muon path, thin lines pointing to the channels mark the path of the Cherenkov 
photons as given by the fit to the measured times. The areas of the circles are proportional to the measured 
amplitudes. Top right: Hit times versus vertical channel positions. Bottom right: The allowed 9/<p regions 
(see text). The fake probability of this event is smaller than 1%. 




Figure 7: a) - an unambiguous 14-hit neutrino candidate; b) - an ambiguous event reconstructed as a 
neutrino event (dashed line) but with a second solution above the horizon (solid line) . This event is assigned 
to the sample of downward going muons. 
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In the mean time, 70 days from NT-96 have been analyzed, and 12 neutrino candidates have been 
found. Nine of them have been fully reconstructed, 3 nearly upward vertical tracks hit only 2 
strings and give a clear zenith angle but ambiguities in the azimuth angle - similar to the two 
events from NT-36. Taking into account the degradation of NT-96 due to failed OMs, this is in 
agreement with MC expectations. 

NT- 200 will have an effective area of ~ 1500 m 2 , after all cuts and averaged over a cone of 60 
degrees about the opposite zenith. It will record about one separable neutrino per day. 



4 Search for Magnetic Monopoles 



GUT monopoles may catalyze processes violating baryon number conservation ||. For reasonable 
velocities f3 < 10~ 3 , a cross section a c = 0.17 • <r //3 2 is predicted for monopole-proton interac- 
tions |jT0|| , with a Q being of the order of magnitude typical for strong interactions. The distances 
between sequential proton decays along the monopole track in water can be as short as 10~ 2 - 10 1 
cm. 

In order to search for GUT monopoles, a special trigger was implemented which selects events 
with short-time increase of the counting rate of individual channels, as expected from sequential 
Cherenkov flashes produced by the proton decay products along a monopole track. During standard 
data taking runs, the monopole trigger condition was defined as > 3 hits within a time window of 
500 fisec in any of the channels. 

The data taken from April 16th to November 15th 1993 with NT-36 were used to search for 
monopole candidates. We requested that one channel had counted > 7 hits and the second channel 
looking to its face and located 7.5 m away along the string > 3 hits during the same 500 ^sec. 
This reduces the number of the experimentally observed candidates from 3.5 • 10 7 to zero. From 
the non-observation of monopole candidates we obtain the upper flux limits (90 % CL) shown in 



Fig. 8 together with our earlier results, limits from 1MB [11| and Kamiokande [|12| and with the 
astrophysical Parker limit. Limits of 10 -15 and 2.7 ■ 10~ 16 cm~ 2 s^ 1 have been obtained by the 
MACRO [14|] and Baksan [13] experiments, respectively. With the whole statistics taken since 1993, 
the minimal detectable flux can be lowered by more than an order of magnitude. 




Figure 8: Upper limits (90 % CL) on the flux of magnetic GUT monopoles as a function of their velocity 
/3, for different catalysis cross sections <r . 
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Apart from GUT monopoles catalyzing baryon decay, relativistic monopoles may be detected as 
well. The large magnetic charge of a monopole (g = 68. 5e) results in a giant Cherenkov radiation, 
equal to that of a 14-PeV muon for a monopole with /3 mon ~ 1- The non-stochastical nature of 
the Cherenkov emission by relativistic monopoles (contrary to a 14 PeV muon!) may be used 
to select monopole candidates. Actuallyeven monopoles with velocities below their Cherenkov 
threshold {(3 m(m ~ 0.75) may be detected, namely by Cherenkov radiation of 5-electrons ( down to 
Pmon ~ 0.6). The effective area of NT-200 for monopoles with (3 « 1 is estimated as 20,000 m 2 . 

5 Limnology 

Apart from its function as a neutrino telescope, the Baikal detector can be used to monitor wa- 
ter parameters. The array permanently records photomultiplier counting rates, and periodically 
parameters like optical transmission at various wavelengths, temperature, conductivity, pressure 
(CDT sondes), and speed of sound. These measurements form a unique data set which can be 
related to CDT measurements at other locations in order to build a comprehensive picture of water 
exchange processes in the lake. 

For reasons of illustration, we show below counting rate variations at various time scales recorded 
with NT-36 in 1993/94. Counting rates of single PMTs as well as the "local trigger rates" (coinci- 
dence rates of the 2 PMTs of a pair) are dominated by water luminescence. 
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Figure 9: Average counting rate of OMs vs. time, compared to bacteria concentration at surface. 

Fig. 9 gives the counting rate over 2 years and compares it to the bacteria concentration measured 
at a distance of 50 km to the NT-200 site, at 10 m depth below surface. In August/September 
we observe an increase of the luminosity to extremely high levels. The changes of the local trigger 
rate are not reflected in the muon trigger rate, since the muon trigger is essentially dominated 
by atmospheric muons, with negligible contribution by random hits (water luminescence or dark 
noise). This is demonstrated in fig. 10 on a shorter time scale, for a time interval of marked changes 
of the local trigger rate following a strong storm at August 3rd, 1993, which had washed a lot of 
water from a nearby river to the lake. 

Fig. 11 shows a short period of about 8 hours when the counting rates sequentially increased, starting 
with the highest OMs and ending with the lowest. From the time shift of the 3 curves a vertical 
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current of 2.3 cm sec 1 is deduced. This is noticeable since the vertical speed of water currents in 
spring, when the vertical water renewal is considered to be most intensively, is only 0.2-0.3 cm sec -1 ! 
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Figure 10: a) Local trigger rates for channel 7 (downward facing) and channel 8 (upward facing) for August 
lst-9th, 1993. The counting rates are averaged over 30min. b) Muon trigger rates (condition 4/1) for 
channel 7 and 8. Counting rates are averaged over 50min. 
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Figure 11: Counting rate of three OMs along one string during an 8 hour interval at Sept. 24, 1993. 



6 Conclusions and Outlook 

The Baikal detector is well understood, and hrst atmospheric neutrinos have been identified. Their 
total number agrees well with Monte Carlo expectations. Furthermore, limits on the flux of GUT 
magnetic monopoles have been derived. The detector has the potential to study so different ques- 
tions like the search for neutrinos from WIMP annihilation in the center of the Earth on the one 
hand, and monitoring of deep water ventilation processes on the other hand. 
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The Baikal site is competitive to Oceans: The stronger absorption may be taken into account by 
a somewhat denser spacing which, on the other hand, might be a natural approach if one focuses 
to lower thresholds than in Oceans. The external noise is of similar magnitude like in Oceans, but 
with strong seasonal variations. The smaller depth has been shown to be no serious drawback, since 
with appropriate methods neutrinos can be separated effectively. The most remarkable advantage 
of the site is the ice cover which allows reliable and inexpensive deployment and retrieval - still an 
practically unsolved problem in the case of Ocean projects. 

After 144 Optical Modules have been deployed in March/ April 1997, the NT-200 detector will be 
completed in April 1998. 

In the following years, NT-200 will be operated as a neutrino telescope with an effective area 
between 1000 and 5000 m 2 typically, depending on the energy. This corresponds, after all cuts, to 
about 1 atmospheric neutrino per day. Presumably still too small to detect neutrinos from AGN 
and other extraterrestrial sources, NT-200 can be used to push the flux limits for neutrinos from 
WIMP annihilation and for magnetic monopoles. It will also be a unique environmental laboratory 
to study water processes in Lake Baikal. 

Apart from its own value, NT-200 is regarded to be a prototype for a telescope 20-50 times larger. 
With 2000 OMs, a threshold of 10-20 GeV and an effective area of 50,000 to 100,000 m 2 , this 
telescope would have a realistic detection potential for extraterrestrial sources of high energy neu- 
trinos. With its comparatively low threshold, it would fill a gap between underground detectors 
and planned high threshold detectors of cube kilometer size. 
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